Exploiting breakdown of the similarity relation for diffuse light transport: 
simultaneous retrieval of scattering anisotropy and diffusion constant 
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As manifested in the similarity relation of diffuse light transport, it is difficult to assess single 
scattering characteristics from multiply scattered light. We take advantage of the limited validity 
of the diffusion approximation of light transport and demonstrate, experimentally and numerically, 
that even deep into the multiple scattering regime, time-resolved detection of transmitted light 
allows simultaneous assessment of both single scattering anisotropy and scattering mean free path, 
and therefore also macroscopic parameters like the diffusion constant and the transport mean free 
path. This is achieved via careful assessment of early light and matching against Monte Carlo 
simulations of radiative transfer. 



Media having appearance and optical properties to a 
large extent determined by multiple scattering of light are 
truly abundant. This is reflected by an enormous mul- 
tidisciplinary interest in the theory of light propagation 
in turbid media and experimental methods for accurate 
assessment of optical properties, structure and chemical 
composition of such systems. Applications range from 
astronomy atmospherical science Q and climate re- 
search [1] to medicine Q, materials science @ and fun- 
damental photonics [(| . Despite being a field with a very 
long history, the complexity of radiative transfer has ren- 
dered approximative treatment prevalent, and there are 
still strong developments regarding accurate theoretical 
descriptions 0, Q ■ Generally, the similarity relation for 
diffusive light transport [j^liol] tells us that it is difficult to 
infer single scattering anisotropy (g) and scattering mean 
free path (£ s ) from multiply scattered light - macroscopic 
transport is essentially governed by the transport mean 
free path, it = 1 S /{1 — g). Measurement of microscopic 
properties are therefore normally made by studying light 
transport close to a source or transport through small 
samples (see, e.g.,[H[lH and references therein). 

Here, we take direct advantage of the limited valid- 
ity of the diffusion approximation of light transport and 
show, numerically and experimentally, that time-resolved 
experiments allow assessment of single scattering char- 
acteristics even when signals are deep into the multi- 
ple scattering regime. Despite that time-resolved diffuse 
spectroscopy is extensively used to assess optical proper- 
ties .structure and chemical composition of turbid media 
[12H24| , as well as for fundamental investigations of light 
transport [25-3Cj, these opportunities has remained un- 
explored. Normally, the technique is considered to op- 
crating in the diffusion regime, and a major theoretical 
and experimental issue has instead been to understand to 



what extent diffusion theory is applicable [3 lM36f | . Even 
in cases where diffusion breaks down and refined mod- 
eling has been used, evaluation has been limited to pa- 
rameters accessible in the diffusion regime (351 - 140} . Since 
single scattering properties often are of great diagnostic 
value in optical characterization, an extension of capabil- 
ities to include also these is important. 

In radiative transfer, transport is modeled using an ab- 
sorption coefficient /i a , a scattering coefficient fi s and a 
scattering phase function. The famous transport equa- 
tion follows from energy conservation 0, 0] • The deriva- 
tion is phenomenological, reducing the complex wave 
propagation into a Poissonian random walk with expo- 
nentially distributed steps between scattering events, and 
with a correlation in directionality given by the phase 
function (cf. @ for a treatment of the relation between 
Maxwell's equations and radiative transfer). The aver- 
age step between scattering events is the scattering mean 
free path, £ s = 1/ fi s . On a macroscopic scale (i.e., after 
many scattering events), the transport equation reduces 
to a diffusion equation. Assuming homogenous random 
distribution and orientation of scatterers, the diffusion 
constant is 
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where g is the average cosine of the scattering angle, and 
ve the energy velocity 0, |42j]. Here we have reached 
the similarity relation mentioned earlier - as long as the 
quantity l t = i s /(l — g) is conserved, macroscopic trans- 
port is similar. It is thus possible to map a diffusive 
process onto a random walk of exponentially distributed, 
isotropic and independent steps of average length £ t . 
The similarity has also been exploited to, e.g., relax the 
knowledge required on g (and phase function) when using 
Monte Carlo simulations (MC) to evaluate time-resolved 
diffuse spectroscopy [35l - l39l . |43| . Having this said, the 
possibility to resolve the ambiguity related to the simi- 
larity relation by looking at early time-of-flights (TOF) 
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was discussed alread y in the classic paper of Patterson, 
Chance and Wilson [1J|. In addition, there are several 
contributions on how q impacts the validity of the dif- 
fusion approximation HH H(J l4^ - |45| , illustrating the 
intimate link between break in similarity and diffusion 
breakdown. Interestingly, it is widely assumed that diffu- 
sion theory is most accurate for g — 0, while recent elab- 
oration shows that the fastest convergence to a Gaussian 
spatial distribution, i.e. diffusion, occurs for g w 0.7 for 
the particular case of the Henyey-Greenstein phase func- 
tion (45[. In any case, the breakdown of diffusion theory 
has most often been looked upon as a complication rather 
than an opportunity. In fact, it often recommended to 
exclude early times during analysis, partly due to the 
breakdown of diffusion theory, partly due to experimen- 
tal uncertainties in timing [18| . While this may be wise 
under certain circumstances, our work aims at demon- 
strating that it is possible to go beyond that perspective. 

We start with a theoretical view on the possibility to 
infer g from TOF distributions. Although analytical so- 
lutions to the radiative transfer equation are becoming 
available @, 5H, [47j , we will rely on Monte Carlo simula- 
tions (MC). MC is a straight-forward and well established 
method to solve the transport equation, capable of han- 
dling arbitrary boundaries, phase functions and source 
terms. When adhering to the new trend of using graphi- 
cal processing units (GPUs) to parallelize computations 
48], simulations are done rapidly. A simulation involv- 
ing injection of 10 8 random walkers injected into a slab 
of thickness L — 8£t, as done here, takes on the order 
of one second on a modern GPU. Our simulations were 
conducted on a GPU using the code in (48j |. but modi- 
fied for a slab geometry. We set absorption to zero and, 
as customary, used the Henyey-Greenstein phase func- 
tion. We start by considering the case of an internal 
refractive index of 71^ = 1.5 and an external of n e —l. The 
influence of g on TOF distributions for total transmis- 
sion is shown in Fig. [TJ The dependence on g is evident, 
but it is also clear that dynamics fall under similarity 
at longer times. Systems with strong forward scattering 
(.9 > 0.7) exhibit similar dynamics also for early times, 
meaning that the similarity relation breaks down when 
moving towards isotropic scattering (this in agreement 
with earlier works, cf. [43, HU). In this case L/£ t — 8, 
so diffusion theory can be expected to rather well de- 
scribe the long-time behavior [34|. The time-integrated 
transmission was very close to 27.0% for all values of 
g , in good a g reement with the prediction of diffusion 
theory [H, Hfl, T D = (£ t + z e )/(L + 2z e ) = 26.6%. 
Similarly, in terms of the ballistic time tt, the decay 
time constant r varied from 5.70 th (<7 = 0.9) to 5.78 tt 
(9 = 0.0), which is in fair agreement with diffusion the- 
ory, t d = (L + 2z e ) 2 /(7T 2 D) = 6.26 t b . 

Let us now turn to our experimental investigation. 
We have studied the time-resolved transmission of light 
(A = 810 nm) through a 252 urn thick turbid slab consist- 
ing of 1 vol% nano-particles embedded in a polymer with 
n=1.5. The slab is located between glass slides, meaning 
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FIG. 1. Simulated time-resolved total transmission for a tur- 
bid slab with L/£ t = 8 for g = 0.0, 0.1, 0.9 (rii=1.5, n e =l). 
Clearly, even deep into the multiple scattering regime, early 
dynamics can exhibit a significant dependence on g. This sim- 
ilarity breaking is clear evidence of diffusion breakdown. That 
a larger g (more forward scattering) reduce that amount of 
early light is related to that the unscattered intensity decays 
with depth z as exp(— z/l s ), and that having g approach one 
while keeping l t fixed is accompanied by a reduction of l a . 



that we are working with index-matched boundary condi- 
tions (n e =iii=l.5). The used scatterers (coated titania- 
particles, Huntsman Tioxide R-XL, USA) have a nominal 
diameter of ^280 nm and an average refractive index of 
about 2.4. Mie calculations show that the single scatter- 
ing anisotropy g is around 0.6. The optical setup used for 
investigations of this sample is illustrated and described 
in Fig. Note that, contrary to the total transmission 
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FIG. 2. Ultrashort pulses at A = 810 nm from a titanium- 
sapphire laser are injected in a small spot on the sample, and 
the central part (30 (ixn spot, monitored by CCD) of the trans- 
mission is collected and resolved in time using an ultrafast 
gating technique. Femtosecond time-resolution is achieved by 
overlapping the transmitted pulse with a fs reference pulse 
on a BBO crystal and monitoring sum-frequency generation 
at 532 nm (a variable delay line gives access to the full TOF 
distribution). 

configuration considered above, only central light is col- 
lected here. The approach described in this work is, in 
fact, very general and applies well to different configura- 
tions (this is further elaborated below). 

Evaluation of obtained experimental data against a 
database of MC simulations for different levels of it and 
g (using only centrally transmitted walkers, in order to 
match experiment), with resolution 0.5 (xm and 0.1 re- 
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spectively, reveals an optimal match for g = 0.6 and 
it = 39 um (i.e., i s — 15.6 \im). Fig. [3] shows the 
experimentally recored TOF distribution along with sim- 
ulations of different levels of g for £ t = 39 pm, as well 
as the result of least squares matching. Fit merits in- 
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FIG. 3. Evaluation of experimental data based on MC. The 
left graph shows simulations for It = 39 urn and different 
levels of g (black lines, g — 0.1, . . . , 0.9) and compares it with 
the experimental data. Experimental data is normalized to 
overlap with MC at late times where difference due to g are 
negligible. The right graph shows obtained fit merits for the 
three best-matching g, revealing that the best least-squares 
fit is obtained for g = 0.6 and £ t = 39 pm. The fitted It is 
in good agreement with measurement on thicker slabs of the 
same material, and the fitted g is in good agreement with Mie 
calculations. 

dicate that the precision is better than ±0.1 for g, and 
±1 pm for if The marked difference in ^-dependence 



between Fig. [T]and Fig. [3] is related to that the former 
considers total transmission, while the latter concerns 
detection of light from a small area around the optical 
axis. Throwing away a large part of the diffuse transmis- 
sion naturally result in a different balance between the 
ballistic and quasi-ballistic component and the multiply 
scattered part. Note also that L/i t is about 6.5, so es- 
tablished rules of thumb 34] implies that diffusion theory 
should fail to describe even late time behavior. However, 
the late-time decay of the total transmission obtained 
from MC (t=3.6 ±0.1 ps) is in good agreement with the 
predictions of diffusion theory (td=3.60 ps). The reason 
behind this agreement is that the validity of diffusion the- 
ory decreases with the refractive index contrast [34| , and 
we are now studying the forgiving case rij = n e . 

In conclusion, we have shown that single scattering 
characteristics such as g and £ s are accessible also after 
long propagation distances. This is potentially a useful 
development, and we hope that our findings will stimu- 
late further investigations of how to proceed beyond the 
similarity relation. The accuracy and precision of the ap- 
proach depends on the transport regime (e.g. L/i t ) and 
the details of angular collection and conversion (cf. 32] ) , 
being an important matter to be addressed in future re- 
search. Finally, it should be noted that measurements do 
not require ultrafast setups per se, as required temporal 
stability and resolution depends on the scale of the sys- 
tem. We therefore anticipate widespread applicability of 
the approach, from characterization of turbid liquids and 
porous media to biological tissues. 
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